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Abstract 

In the context of the minimal supersymmetric seesaw model, we study the implications 
of the current neutrino data for thermal leptogenesis, (3(3q v decay, and leptonic flavour- and 
CP-violating low-energy observables. We express the heavy singlet-neutrino Dirac Yukawa 
couplings {Y u )ij and Majorana masses Mjv 4 in terms of the light-neutrino observables and 
an auxiliary hermitian matrix H, which enables us to scan systematically over the allowed 
parameter space. If the lightest heavy neutrino Ni decays induce the baryon asymmetry, 
there are correlations between the Mn 1 , the lightest active neutrino mass and the primordial 
lepton asymmetry e\ on the one hand, and the /3/3o u decay parameter m ee on the other hand. 
However, leptogenesis is insensitive to the neutrino oscillation phase. We find lower bounds 
Mn 1 10 10 GeV for the normal light-neutrino mass hierarchy, and <: 10 11 GeV for 
the inverted mass hierarchy, respectively, indicating a potentially serious conflict with the 
gravitino problem. Depending on M^, we find upper (upper and lower bounds) on the 
lightest active neutrino mass for the normal (inverted) mass hierarchy, and a lower bound 
on m ee even for the normal mass ordering. The low-energy lepton-flavour- and CP-violating 
observables induced by renormalization are almost independent of leptogenesis. The electron 
electric dipole moment may be close to the present bound, reaching d e ~ 10~( 27-28 ) e cm in 
our numerical examples, while may reach d^ ~ 10~ 25 e cm. 



CERN-TH/2002-133 
June 2002 



1 Introduction 



The only convincing experimental evidence for physics beyond the Standard Model, so far, is 
provided by neutrino oscillations fl], |2|, which are generally interpreted as evidence for non- 
zero neutrino masses. The smallness of neutrino masses is explained naturally via the seesaw 
mechanism || with heavy singlet Majorana neutrinos. Leptogenesis scenarios envisage that 
CP violation in their out-of-equilibrium decays may induce a non-zero lepton asymmetry, 
which may be converted into the observed baryon asymmetry of the Universe via electroweak 
sphaleron processes [|J. If, in addition, there is low-energy supersymmetry as motivated by 
the hierarchy problem, which is exemplified by the vastly different mass scales involved in 
the seesaw mechanism, low-energy processes violating charged-lepton flavour and CP may 
be observable, induced by the neutrino Dirac Yukawa interactions via renormalization below 
the heavy-singlet mass scale || |J. In this scheme, the only source of many observables 
- neutrino oscillations, PPq u decay, the baryon asymmetry of the Universe, and flavour- 
violating decays of charged leptons and their electric dipole moments (EDMs) - are the 
Dirac Yukawa couplings (Yv)ij and the Majorana masses of the three heavy singlet 
neutrinos. 

Today, most of our knowledge of Y u and comes from the light-neutrino mass and 
mixing parameters measured in oscillation experiments, with an additional constraint from 
searches for f3f3o u decays 0]. The oscillation data are converging towards unique solutions 
for each of the atmospheric and solar neutrino anomalies ||, which represents a major 
breakthrough in neutrino physics. On the other hand, the determination of the baryon 
asymmetry of the Universe has significantly improved in recent years ||, and will improve 
still further with further astrophysical and cosmological observations, such as those of the 
MAP and Planck satellites. In addition, one expects significant improvements in the future 
experiments searching for j3/3 Qu decay fL0| , lepton-flavour- violating (LFV) processes |ll], [12], 



13| and electric dipole moments (EDMs) |14], [U| [OJ. These prospects motivate us to perform 



a comprehensive study of these observables in the minimal seesaw model. 

Any study of leptogenesis, neutrino masses, LFV processes and the EDMs of charged 
leptons faces the generic difficulty of relating the experimental information on light neutrino 
masses and mixings with other observables. If one takes a top-down approach and fixes Y u 
and Mjy. by some theoretical or phenomenological argument [[16], such as GUT relations, U(l) 
or non-Abelian flavour models, phenomenological textures, democratic principles, arguments 
of minimal fine-tuning, etc., one can study the pattern of typical predictions of the model 
considered, but cannot perform a comprehensive phenomenological study of the interesting 



1 



observables. Even correct numerical consistency with light-neutrino data is a difficult task 
in this approach, since it may involve fine-tunings and must be checked a posteriori. 

These problems can be solved in a bottom-up approach [|17]] to neutrino observables. 
Parametrizing (Y v )ij and in terms of the light neutrino mass matrix M. v and an auxil- 
iary Hermitian matrix H, as in (TTS|, compatibility with the light-neutrino data is automatic, 
because Ai u is an input. In addition, since the Hermitian matrix H has a physical inter- 
pretation as H = YjDY u , where D is a real diagonal matrix, one has also control over the 
renormalization-induced LFV decays. For every (Y u )ij and M^. generated in this way, one 
can therefore calculate exactly the weighted light neutrino mass m ee measured in f3(3Q U de- 
cay [19], rates for LFV decays |6], 20, 21], and EDMs of charged leptons |22| ^3fl . Moreover, 



one can also calculate consistently the leptogenesis CP asymmetries e; t [24], ^5], and, assum- 



ing the standard thermal leptogenesis scenario, also the washout parameters |26 



Hence one can study the correlations between all these parameters and their dependences 
on the light and heavy neutrino masses. To our knowledge, there has so far been no study 
of leptogenesis in which all the e^, Hi, m v . and M^. are treated simultaneously as dynamical 
variables determined in consistency with the oscillation data. 

In this paper we perform such a comprehensive phenomenological study of three classes 
of leptonic observables in the minimal supersymmetric seesaw model: the effective light- 
neutrino parameters in the baryon asymmetry of the Universe generated via thermal 
leptogenesis, and the renormalization induced LFV processes and EDMs of charged lep- 
tons. Assuming the large-mixing-angle (LMA) solution to the solar neutrino anomaly, we 
parametrize Y v and M^r. in terms of Ai u and H. In our previous paper [18| we considered 
only the case H = YjDY u where Da = l^Mciir/M^). Here we also consider a different 
form for D |~7|, which yields maximal EDMs for the charged leptons. We assume the stan- 
dard thermal leptogenesis scenario in which the observed baryon asymmetry of the Universe 
is generated only by the decays of the lightest singlet neutrino iVifl. The wash-out parame- 
ter K\ is calculated by solving numerically the analytical Boltzmann equations of [29], and 
correcting by a constant factor in order to be consistent with the exact results in J28|. 

We scan randomly over the input parameters: the lightest neutrino mass m Ul (or m U3 for 
the inverted hierarchy of light-neutrino masses), the two Majorana phases in Ai u , and the 
entries of the parameter matrix H. We require the Yukawa couplings Y v to remain perturba- 
tive until the GUT scale, we require the generated baryon asymmetry to be consistent with 



1 This assumption is not valid in the case of highly-degenerate heavy neutrinos, but still holds well in the 
case of moderate degeneracy [25 . As the are output parameters in our approach, it is not suitable for 
systematic studies of very degenerate heavy neutrino phenomena. 
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observation, and we also require consistency with all the bounds on LFV processes. 

We find interesting correlations between the neutrino observables and leptogenesis pa- 
rameters, whilst the LFV processes and EDMs are almost independent of the leptogenesis 
constraints. The influence of the neutrino oscillation CP phase 5 on leptogenesis is negli- 
gible: the existence of a baryon asymmetry does not require it to be non-vanishing. The 
experimental bound Yg ^ 3 x 10~ n on the baryon-to-entropy density ratio Yb implies the 
lower bounds Mjv x ^ 10 10 GeV and ^ 10 11 GeV for the normally and inversely ordered 
light neutrino masses, respectively 0. These bounds put the findings of [|30) on rigorous 
numerical ground, and indicate a serious potential conflict with the gravitino problem in 



supergravity models |[1 1 . Leptogenesis also implies non-trivial bounds on the mass of light- 
est light neutrino. For the normal mass ordering, there is an -dependent upper bound 
on m vi , whilst for the inverted hierarchy there are both upper and lower bounds on m U3 . 
Successful leptogenesis with m Ul ^ 0.1 eV could be allowed for ^ 10 12 GeV. There is 
also an -dependent lower bound on m ee for normally-ordered light neutrinos, implying 
its possible measurability in future experiments. On the other hand, m ee has a preferred 
value determined by Am 2 tm even in the case of the inverted mass hierarchy. It tends to be 
below O(10 _1 ) eV, making improbable the discovery of (3(3q v decay in current experiments. 

The rates of LFV processes and EDMs depend also on the soft supersymmetry-breaking 
parameters, which we fix by choosing one of the post-LEP benchmark points . We find 
that Br(r — > fi{e)^f) and Br(fi — > ej) may easily saturate their present lower bounds, and 
that the EDMs of electron and muon may reach d e ~ icr( 27 ~ 28 ) e cm and ~ 1CT 25 e cm, 
respectively, in our random samples. We stress in particular that the electron EDM d e may 
be less than an order of magnitude from the present experimental bound d e ^ 1.6 x 10 -27 
e cm |33|] and offers a sensitive probe of the supersymmetric seesaw model. We find also 
some correlation between LFV r decay rates and the lower bound on . 

The paper is organized as follows. In Section 2 we classify the observables and discuss our 
parametrization. Section 3 contains phenomenological studies, and finally our conclusions 
are presented in Section 4. 

2 Parametrization of Neutrino Observables 



2 These conclusions on neutrino parameters and leptogenesis are valid also in non-supersymmetric models. 
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2.1 Observables and Physical Parameters 

The leptonic superpotential of the minimal supersymmetric model which implements the 
seesaw mechanism is 

W = N^{Y v ) l3 L 3 H 2 - Ei(Y e )ijLjHi + ^N c l {M N ) lJ N'r + ^H 2 H X , (1) 

where the indices i,j run over three generations and is the heavy singlet-neutrino mass 
matrix. One can always work in a basis where the charged leptons and the heavy neutrinos 
both have real and diagonal mass matrices: 

(Y e ) l3 = YgSij , (M N )ij = M Ni 5ij . (2) 

The matrix Y v contains six physical phases and can be parametrised as 

(Yv)ij — Z ik Y®Xlj, (3) 

where X is the analogue of the quark CKM matrix in the lepton sector and has only one 
physical phase, and Z = PiZP 2 , where Z is a CKM- type matrix with three real mixing 
angles and one physical phase, and Pi 2 = diag(e 1,3 J e 2A , 1). This implies that we have 15 
physical parameters in the Yukawa coupling Y u , which together with the 3 unknown heavy 
masses Mat make a total of 18 parameters in the minimal seesaw model []. 

These 18 unknown neutrino parameters give rise to three classes of physical observables, 
as presented diagrammatically in Fig.|l|: 

(i) Low-energy effective neutrino masses arising from the seesaw mechanism: 

M u = Yj {M N )~ l Y u v 2 sin 2 (3. (4) 

The effective light-neutrino mass matrix M. v is symmetric and can be diagonalized by 
a unitary matrix U as follows: 

U T M V U = M°. (5) 

By a field redefinition, one can rewrite U = VPo, where Pq = diag(e*^ 1 , e 1 ^ 2 , 1) and V is 
the MNS matrix. Therefore, all the low-energy neutrino observables such as neutrino 
oscillations, (3(3q v decay, etc., depend on the 9 effective parameters in Ai u , which are 
functions of all the 18 parameters in ([!]). Whilst neutrino oscillations measure the 
3 The parameter counting is identical in models with and without supersymmctry. 
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Seesaw mechanism 

M v 

9 effective parameters 




Leptogenesis 
Y„Yt , M Ni 
9+3 parameters 



Renormalization 

YjLY, , M Ni 
13+3 parameters 



Figure 1: Roadmap for the physical observables derived from Y v and iVj. 



mass-squared differences of neutrinos and their mixing angles, PPq v decay measures 
one particular combination of their masses and mixing matrix elements, 



(6) 



which involves also the Majorana phases. The NMS mixing phase 5 can in principle be 
measured in neutrino oscillations with neutrino factory experiments, but measurements 



of the Majorana phases are less straightforward |19| |34 



(ii) The idea of baryogenesis via leptogenesis is first to produce a lepton asymmetry, whose 
ratio to the entropy density we denote by Yl, via the out-of-equilibrium decays of the 
heavy neutrinos iVj. This asymmetry is converted into the baryon asymmetry of the 
universe via (B + L)-violating electroweak sphaleron processes: 



Yn 



C 



c 



-Yr 



(7) 
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where C = 8/23 in the minimal supersymmetric extension of the Standard Model 
(MSSM). The generated lepton asymmetry depends on the initial neutrino-plus-sneutrino 
density, the CP asymmetries in neutrino decays, and on the washout effects. The CP 
asymmetry produced in out-of-equilibrium decays of iVj is given by |J§] 



Im 




[Y U Y/ 


) 




il - 




1 V 


2 



v^[Log(l + 1/xi) + 2 _ ] , (8) 



where x\ = (Mjv,/MjyJ 2 . It is clear from @ that the generated asymmetry depends 
only on the 9 parameters (including 3 phases) in 

Y v Yl = P{T{Y») 2 Z T P X (9) 

and on the heavy neutrino masses. 

In the case of non-degenerate heavy neutrinos, the lepton asymmetry is, to a good 
approximation, generated only via the decays of the lightest heavy neutrino N\ (and 
the corresponding sneutrinos), because the very rapid washout processes mediated by 
N\ erase whatever asymmetry had been produced previously in jV 2) 3 decays ]26|, j2~5| . 
In this case one has a simple relation 

y L = y J 5?(o)ei«i, (io) 

where Y^ q (0) is the initial thermal equilibrium density of the lightest neutrinos and k% is 
the washout parameter. The factors multiplying e\, in ( |10"D depend on the cosmological 
scenario. In the case of the standard leptogenesis scenario with thermally produced 
heavy neutrinos one has in thermal equilibrium the initial condition Y^(0) ~ 1/(2^*), 
where g* ~ 230 is the number of effective degrees of freedom in the MSSM. The washout 
parameter k± can be precisely calculated by solving the set of Boltzmann equation for 
Yl and for the neutrino density Y^, 



dY Nl 
dz 



+ 7iVi(^) = 



-T L + 'Yl(z)Y l + S l (z)=0, (11) 
dz 

where Y L , Y Nl and the factors 7^(2), 1l{ z ) and Sjj(z) depend on the temperature, 
which is parametrized by z = M 1 /T. In this scenario, 7^(2), Sl(z) and 7l(-z) are 
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worked out in great detail [^6], 28] in the Standard Model, for which analytical approx- 

In general, the 



imations exist P9| , and also in the supersymmetric framework 
generated asymmetry depends on four parameters: the CP asymmetry e±, the heavy 
neutrino mass Mfj x , and the lightest light neutrino mass, which determines the overall 
light-neutrino mass scale. The effective mass parameter 



mi 



v 2 sin 2 j3 



M 



(12) 



is also used in discussions of leptogenesis: it depends on the 10 parameters introduced 
above. Fixing the baryon asymmetry to agree with observation implies a relation 
between e\ and K\ via (]T0|) . Therefore, only the relations between e\ and the masses 
Mtvj and m Ul (or m U3 for the inverted mass hierarchy) are physically relevant. 

(iii) Renormalization of soft supersymmetry-breaking parameters due to the presence of 
Y v above the heavy-neutrino decoupling scales induces low-energy processes such as 
the charged-lepton decays // — > ej, \i — > eee, r — > Ij, r — > 3/, and (in the presence 
of CP violation) EDMs for the electron and muon. In a leading-logarithmic approx- 
imation, renormalization modifies the left-slepton mass matrix and trilinear soft 
supersymmetry-breaking A e terms according to || 



(5m 

which are proportional to 



1 



(3ml + Al)^LY) 



8tt 2 

l -A Q Y ei {Y^LY) ij: 



Stt 



Y^LY 



xy d p 2 z t lz*p;y d x\ 



where L is a diagonal matrix 



\n{M GUT /M Ni )8,. 



(13) 



(14) 



(15) 



We note that the CP-violating observables at low energies depend on the leptogenesis 



phase in Z. The expression (|T4| ) contains 16 neutrino parameters altogether, but in com- 
pletely different combinations from the seesaw mass matrix (01). Only the two Majorana 
phases in Pi cancel out in (0). If the heavy neutrinos are exactly degenerate in mass, 
all the renormalization-induced observables are proportional to YJY V = X{Y V D ) 2 X\ 
but this is not a good approximation, in general. 
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The CP- violating observables in LFV processes all depend on a single CP invariant, 
which is related to H = YjLY u by J = Imifi 2 i/23-f^3i [051 . This influences slepton 



physics at colliders and also determines the T-odd asymmetry in \x — > 3e J37|. The 
dominant contribution to the lepton EDMs arises from threshold corrections to the 
trilinear coupling A e due to the non- degeneracy of heavy neutrino masses. Diagonal 



phases in A e are proportional to [|3|, |18 



ImlXj, X k ]u log M N JM Nj ^ 0, (16) 

where (-Xfc)y = {Y* (M^)) ki(Y u (MN k )) kj ■ This depends non-trivially on the CP-violating 
phases, including the two Majorana phases in M. v and two phases in H that are irrel- 
evant for LFV. 

2.2 Phenomenological Parametrization of Y v and M^. 

None of the three classes of observables discussed in the previous subsection, by itself, al- 
lows one to measure all the parameters in the neutrino superpotential (|l|). Moreover, at the 
moment the only parameters known experimentally are the 2 effective light neutrino mass 
differences and 2 mixing angles, and the baryon asymmetry of the Universe, whose interpre- 
tation requires some cosmological inputs. A central issue for comprehensive studies of the 
neutrino sector is how to parametrize Y v and M^r. in such a way that the effective neutrino 
parameters measured in the oscillation experiments are incorporated automatically, and all 
other phenomenological constraints are also satisfied. 

One may attempt to fix Y v and Mjy. in ([!]) using model predictions or some other principle. 
In that case, however, satisfying the measured effective neutrino masses and mixings in (H) is 
a non-trivial task. One can study the patterns of typical predictions for LFV processes, EDMs 
and the baryon asymmetry in any specific model, but not make comprehensive numerical 
studies which cover all the allowed parameter space. 

However, as discussed in [O , in the supersymmetric seesaw model the low-energy degrees 
of freedom may in principle be used to reconstruct the high-energy neutrino parameters. 



In [[Tj| we presented a parametrization of Y v and Mj^ i in terms of the light-neutrino mass 



matrix M. u and an auxiliary Hermitian matrix H, 

H = YjDY u , (17) 



where the diagonal matrix D was chosen in [18| to be Dij = \n{MGUT/M^.)8ij^ moti- 



vated by (|14D . In this case, the parameter matrix H is directly related to the solutions 
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of the renormalization-group equations (RGEs) for the soft supersymmetry-breaking slepton 
masses, according to (|IBD, and allows us to control the rates for LFV processes. Conversely, 
if any LFV process is observed, one can use its value to parametrize the heavy neutrino 
masses and couplings. 

In principle the diagonal matrix D can be an arbitrary real matrix. In order to study 
the maximal range for the charged lepton EDMs in the supersymmetric seesaw model, we 
also study in this paper the parametrization with D = 1, which is the case discussed in JL7). 
According to (16), the EDMs are not proportional to H in either case. However, the choice 



D = 1 departs from (16) more than the choice D = L. As some of the entries of H must 
be suppressed in order to satisfy the stringent experimental upper bound on \i — > cy, one 
expects the EDMs to get somewhat larger values if D — 1, and this is indeed the case. 
We now present the details of the parametrization. Starting with ([17]) and using the 



parametrization 21 



v sin p 



one can recast (IlTT) into a form 



H' = R^M N R\ (19) 



where is a diagonal matrix 



and 



(M N )u — D U M N ., (20) 



H'=^M° U^HU^M® w 2 sin 2 /?. (21) 

If (plj) can be solved, i.e., if the matrix H' can be diagonalized with an orthogonal matrix 
R f , then one can solve the heavy neutrino masses from ( |20"D and calculate the neutrino Dirac 
Yukawa couplings from ([18]). Schematically, 



(Mu , H) — > (M v , M N , R') — > (Y„ ,M Ni ), (22) 



where the quantities and R' are calculated in the intermediate step and do not have any 
independent physical meaning. Thus, one has converted the 9 low-energy effective neutrino 
parameters and the 9 free parameters in the Hermitian parameter matrix H into the 18 
physical neutrino parameters in Y v and M^ v 
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The advantages of this parametrization have already been mentioned: it allows one to 
control the rates for LFV processes and to scan efficiently over the allowed parameter space 
at the same time. The disadvantage of the parametrization ( p2| ) is that it is not continuous, 
because for some choice of the parameter H there may not exist a matrix R! that diagonalizes 
H' . However, in the case of mult i- dimensional parameter spaces such as the one we study, 
scanning randomly over the allowed parameters is the most powerful tool, and in practice 
this disadvantage does not hinder such a phenomenological study 



3 Phenomenological Analysis 

Using the parametrization developed in the previous section, we now perform a compre- 
hensive phenomenological study of all three types of leptonic observables in the minimal 
super symmetric seesaw model. 

We fix the known light neutrino parameters by Am| 2 = 3 x 10~ 3 eV 2 , Am^ = 4.5 x 10~ 5 
eV 2 , tan 2 #23 = 1 and tan 2 #12 = 0.4, corresponding to the LMA solution for the solar neutrino 
anomaly. Since the experimental constraint on the angle #13 is quite stringent, sin 2 2#i3 ^0.1 
p8| . p9| , our results depend only weakly on its actual value. We fix sin 8 13 = 0.1 and study two 
cases with the limiting values of the neutrino mixing phase 5 = 71/ 2 and 5 = 0. We consider 
both normally and inversely ordered light neutrino masses, since neutrino oscillations do not 
discriminate between them at present^. 

We assume the standard thermal leptogenesis scenario in which the baryon asymmetry 
originates only from the lightest heavy neutrino decays, as described by (|T0l). In our cal- 
culations, we solve (|TT|) numerically using the approximate analytical expressions for the 
thermally-averaged interactions given in ||29|| . We start at T ^> with the initial condi- 



tions Fjvi = Yjsrv an d = 0- As i s appropriate for supersymmetric models, we concentrate 
on the Iow-Ma^ regime, ^ 10 13 GeV, where the approximate solutions for Ki differ from 
the exact ones p8[ by just a constant factor. We correct our output by this factor to be 
consistent with [[28| . These results for the washout parameter Ki were derived in the context 
of non-supersymmetric seesaw models, but are expected to be a good approximation also in 
supersymmetric models, especially for the low and the moderate values of tan/3 that 
we consider in this work. Because the heavy neutrinos thermalize very fast, our results are 
actually independent of the initial conditions in this Iow-Ma^ regime. Here the wash-out pa- 
rameter Ki depends practically only on the effective parameter rh\ (with a small dependence 
on MjVi). Therefore, with high accuracy, the observed baryon asymmetry implies via ( |TUD 
4 Future neutrinoless double beta decay experiments could resolve this ambiguity [lJJ. 
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that ei and rh\ have a one-to-one correspondence. 

In our subsequent analysis, we require the induced baryon asymmetry be in the range || 

3 x KT 11 < Y B < 9 x 10 -11 , (23) 

and study its implications on the light and heavy neutrino masses, the CP asymmetry e±, 
the flflvo decay parameter m ee , the LFV decays and the EDMs of the charged leptons e and 
ft. 



As was shown in [18|, the stringent experimental limit on Br(/2 — > cy) implies the phe- 
nomenological constraints #12 <C 0(1) and H13H32 <C 0(1) on elements of our parameter 
matrix H. Therefore we study two different textures of the matrix H: 




#1 = [ b d I (24) 
and 

/ a d \ 

#2 = 6 0, (25) 
V dt c / 

where a, 6, c are real and d is a complex number. The texture #1 suppresses r — >• and d e 
while r — > /I'y and can be large, and vice versa for H 2 , since these processes are sensitive 
to #13 and #23, respectively. We consider two forms of the matrix H (|i~7|), namely those with 
D = L and D = 1. For the textures #1, #2, the former suppresses \x — > very efficiently 
and all the parameters a,b,c, \d\ can be simultaneously of order unity. However, for the 
D — 1 case one of a, b, c must necessarily be very small in order to keep fi — > e7 below the 
present experimental bound. 

In our model, the rates for the LFV processes and EDMs depend on the soft supersymmetry- 
breaking parameters. In the following numerical calculations we fix them at the GUT scale 



to coincide with one of the post-LEP benchmark points m.1/2 = 300 GeV, mo = 100 
GeV, Aq = —300 GeV, tan/3 = 10 and sign(n) = +1. This choice ensures that all other phe- 
nomenological constraints, such as those on the lightest Higgs boson mass, supersymmetric 
contributions to b —> 57 and — 2, cosmological arguments, etc., are satisfied within the 
necessary accuracy. 

As the input parameters we thus have the lightest neutrino mass m ui (or m V3 for inversely- 
ordered neutrinos), the two low-scale Majorana phases 01,2, and the entries a,b,c,d of the 
parameter matrices #1 and #2. We generate these input parameters randomly in the ranges 
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(1CT 4 — 1) eV for m Vl (or m Vz ), (0 — 2tt) for 0i i2 , and (10~ 7 — 10) for a,b,c, \d\. The distribution 
for each of them is flat on a logarithmic scale. We require that Y v remains perturbative up 
to Mqut and impose the present constraints on all the LFV processes. 

3.1 Normally- Ordered Light Neutrinos 

First we study the implications on leptogenesis parameters, /3j3 u o decay, and LFV decays and 
charged lepton EDMs in the case of normally-ordered light neutrino masses and mixings. 
Apart from the EDMs, the two parametrizations with D = L and D = 1 give practically 
indistinguishable results. Therefore, we present the D = 1 results only for the EDMs and the 
LFV decays, while all the other plots present results for the parametrization with D = L, 
as in 0. 




Figure 2: Scatter plot of the CP-violating asymmetry e\ and the lightest heavy neutrino mass 
Mn 1 for the two extreme choices of the MNS phase 5 = n/2 and 5 = 0, assuming normally- 
ordered light neutrinos and the texture H\. The baryon asymmetry is required to be in the 
range (^). 

We present in Fig. |2] scatter plots for the CP-violating asymmetry ei and the lightest 
heavy neutrino mass Mjv 1( for two extreme values of the MNS phase: 5 = tt/2 and 5 = 0, 
assuming the texture Hi. Here we also introduce a colour code to study the distribution of 
points with different in the following plots. The points within a factor of five from the 
lower bound on Mn x are black, while the points with larger are grey (green). 



12 



Fig. ^ shows immediately that there is no distinction between the plots for 8 = it/2 and 
for S = 0. In special cases, there have been studies whether the observed non-zero baryon 



asymmetry can be related to the NMS phase S |[40||. Our results imply that this is not 



possible in general, and that successful leptogenesis does not require a non-zero value for S. 

We also see immediately in Fig. |2]that there is an M Nl -dependent upper bound [^] on the 
cosmological CP- violating asymmetry e\, and there is a strong lower bound 10 10 GeV 

on the Ni mass. This indicates a potential serious conflict with conventional supersymmetric 
cosmology which requires an upper bound on the reheating temperature of the Universe after 
inflation, derived from avoiding gravitino overproduction |31 . 



In [30], the analytical bound 

3 M Nl 
07T v z sin p 

was found in the limit of very hierarchical heavy neutrinos, <C Mn 2 -C Mn 3 . Our 

results improve this bound, by including the best available numerical results for the washout 
parameter k± and allowing for moderately degenerate heavy neutrinos. We have compared 
the bound ( p6|) with our numerical calculations. It is well satisfied for the low-M^ points in 
Fig. [| but can be violated by a large factor for points with high and high e\. In these 
cases, the heavy neutrinos are not hierarchical in mass and e\ is therefore enhanced. As we 
see in the next figures, for these points the light neutrino masses can also be moderately 
degenerate. 

To study correlations between the leptogenesis parameters, we give in Fig. |3| scatter plots 
of the CP-violating asymmetry e\ and the effective mass parameter fh\ versus the lightest 
neutrino mass m Vl for the texture H\. The shape of the both plots is the same, verifying 
that, for fixed Yb , the parameters e\ and rh\ are not independent parameters in our scenario, 
as discussed above. There is no lower bound on m Ul in this scenario, and the upper bound 
on m Vl depends on , as indicated by the distribution of colours. The allowed band for e\ 
depends only weakly on m Vl for m Vl 

10(-!- 2 ) eV, but the depend ence becomes strong for 
larger m vi . While degenerate light neutrino masses are disfavoured by our results, confirming 
the claims in [^, ^T|, light neutrino masses m ui ~ 0(0.1) eV are still perfectly consistent 
with leptogenesis. Also, notice that for small m Vl the lower bound on fhi is much stronger 



than the limit m\ > m u , derived in J28L 41 



To study the possible implications of the light neutrino masses and leptogenesis for /3/3 0u 
decay, the LFV processes and EDMs in the texture Hi, we give in Fig. § a scatter plot of 
the Mjvi against the /3(3q u parameter m ee , and the muon electric dipole moment versus 
the branching ratio for r — > fi'-f. Surprisingly, there are both lower and upper limits on m ee , 
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Figure 3: Scatter plot of the CP-violating asymmetry €\ and the effective mass parameter fh\ 
versus the lightest neutrino mass m Ul for the normal mass hierarchy and the texture Hi. 



io" 4 





m PP [eV] 




BR(x^hy) 



Figure 4: Scatter plots of M Nl and the (5(3q v parameter m ee , and of the muon electric dipole 
moment d M and iae branching ratio for r — > 7x7 decay, for the normal neutrino mass ordering 
and the texture Hi. 
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which depend on M Nl . For M m < 10 11 GeV we get 10 -3 < m ee < 10 -2 eV. 

There is no correlation between leptogenesis and the renormalization-induced LFV decays 
and EDMs, as seen in Fig. f|. With the chosen soft supersymmetry-breaking parameters, 
Br(r — > fi'y) can attain the present experimental bound in our random sample, and the 
muon EDM may reach ~ 10~ 25 e cm. The former may be observable at the LHC and 
in B-factory experiments, which may reach sensitivities Br{r — > fij) ~ lO 1 - -8-9 ^ [42|, and 
the latter in experiments at the front end of a neutrino factory, which may be able to reach 
d>n ~ 5 x 1CT 26 e cm fll3fl . The texture Hi suppresses Br(r — > ej) and d e below observable 
limits. 




10" ' 10"'" 10"" 1(f lu 10-" 10"° 10"' 10"° 

m ee [eV] BR(x^ey) 



Figure 5: Scatter plots of versus the (5(3q v parameter m ee , and of the electron electric 
dipole moment d e versus the branching ratio for t — > e^y, for the normal neutrino mass 
ordering and the texture H 2 . 



We have performed a similar analysis for the texture H 2 . The behaviour of the leptogenesis 
parameters is almost the same as for the texture H±, already shown in Figs. |2| and |3], so we 
do not present further plots here. The most important differences from Hi can be seen, 
however, in Fig. || where we present scatter plots of M Nl versus the j3j3 0u parameter m ee , 
and of the electron electric dipole moment d e versus the branching ratio for r — > e^y. Whilst 
the lower bound on is the same as in the previous case, the distribution of points 
clearly favours large values of . This is a result of a mismatch between the structure of 
H 2 and the light neutrino mass hierarchy m Ul < m U2 < m Vz . Because Hn, if 13 7^ in H 2 , the 
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Yukawa couplings for the first generation tend to be larger than in the case of the texture 
Hi. As the lightest N\ tends to be related to the lightest light neutrino mass m Vl , the seesaw 
mechanism implies that larger M^ t are usually needed. Fig. |5] indicates that this mismatch 
can be compensated with a tuning of the input parameters, resulting in a relatively small 
number of points at low M^ 1 . 

For the texture H2, Br{r — > ey) and the electron EDM can be large, as seen in Fig. [5], 
whilst Br{r — > yq) and are suppressed below the observable ranges. We find that 
Br(r — ► erf) can be of the same order of magnitude as Br(r — > yy), shown in Fig. f|. 
Importantly, the electron EDM may exceed d e ~ 10~ 28 e cm in our numerical examples. 
This is just one order of magnitude below the present bound d e ^ 1.6 x 1CT 27 e cm J33j. As 
we have not made special attempts to maximize it, it might even reach larger values with 
special values of the soft supersymmetry-breaking parameters. 

We see that there are no black points in the right-hand panel of Fig. [5], corresponding 
to the facts that Brir — > eq) and d e are very much suppressed. This implies a correlation 
between the lower bound on Mfj x and the rates of the FLV processes and EDMs. If Brir — > 
ej) were to be found at the LHC while Brir — > yy) were not, the lower bound on from 
leptogenesis would be above 10 11 GeV. 



3.2 Inversely-Ordered Light Neutrinos 

The inverse ordering of light neutrino masses, m U2 > m Vl > m U3 , is still a viable option, since 
present neutrino oscillation experiments are not sensitive to the sign of Ami. Therefore, it 
is interesting to study whether this case can be discriminated from the normal hierarchy of 
light neutrino masses, using leptogenesis, LFV and CP-violating processes. 

In the inverted light-neutrino mass hierarchy, the correlations between ei and the heavy 
and light neutrino masses are practically the same for both textures Hi^, and it is sufficient to 
present results for just one of them. Since d M is suppressed in Hi by two orders of magnitude 
compared with the normal hierarchy, due to smaller third-generation Yukawa couplings, we 
choose H2 for presentation. 

We present in Fig. ^ scatter plots for the lightest heavy neutrino mass and the 
CP- violating asymmetry ei versus the lightest light neutrino mass m U3 . As expected, the 
lower bound on is higher than in the normally-ordered case, and is *S 10 11 GeV. 
This result follows from the seesaw mechanism, since the Yukawa couplings for the first two 
generations must be larger. Somewhat surprisingly, leptogenesis implies also a lower bound 
on the lightest neutrino mass m V3 . Therefore the light neutrinos tend to be degenerate if 
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1^3 [eV] 



Figure 6: Scatter plots of the lightest heavy neutrino mass and of the CP-violating 

asymmetry e% versus m„ s , for inversely-ordered light neutrino masses and the texture H 2 . 




m ee [eV] 




Figure 7: Scatter plots of versus m ee and of d e versus Br(r — > e^y), for the inversely- 
ordered light neutrino masses and the texture H^- 
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their masses are inversely ordered. The CP-violating asymmetry e% can be larger than in the 
previous case, implying stronger washout effects and somewhat larger light neutrino masses, 
as seen in Fig. |6|. 

In Fig. |7|we give scatter plots of versus m ee and of d e versus Br{r — > erj) for the same 
texture as previously. The sharp lower bound on m ee is the consequence of the inverted mass 
ordering. There is a preferred region for m ee determined by Am^ tm , and relatively few points 
extend above m ee ^ 0(0.1) eV. Therefore, even for the inverted mass hierarchy, observation 
of PPov decay in current experiments is improbable. Again, d e and Br[r — ► cy) reach the 
same values as in the case of normally-ordered neutrinos, and no strong correlation between 
leptogenesis, d e and Br(r — > ej) is present. 



4 Discussion and Conclusions 

In the context of the minimal supersymmetric seesaw model, we have studied relations 
between the light and heavy neutrino masses, thermal leptogenesis, and LFV decays and 
EDMs of charged leptons, scanning over the phenomenologically-allowed parameter space as 
suggested in |Tgf. There are lower bounds M Nl > 10 10 GeV and M Nl k, 10 11 GeV for the 
normal and inverse hierarchies of light neutrino masses, respectively. These bounds are in 
potential conflict with the gravitino problem in supersymmetric cosmology if the gravitino 
mass is below TeV. 

In the thermal leptogenesis scenario, one can avoid these bounds by fine-tuning model 
parameters, namely the Yukawa couplings Y v and the heavy neutrino masses M^. It is 
well known that the CP-violating asymmetry may be enhanced if the heavy neutrinos are 



degenerate in mass |25[ . This may allow one to lower [[25] , |43|1 the lightest heavy neutrino mass 
MjVi , and so avoid the bounds on the reheating temperature [ 31 1 . Degenerate heavy neutrinos 
are also consistent with the LMA solution to the solar neutrino problem | |44j| . Because our 
parametrization gives Y u and as an output, such fine tunings cannot be studied by our 
random scan over the parameter space, and would require a different approach. 

Another way out of the problem would be leptogenesis with non-thermally produced 



heavy neutrinos |41| . In such a case, the reheating temperature of the Universe does not 
limit leptogenesis. However, they do not have the same predictivity, and the implications 
of the observed Y B on neutrino parameters are lost. It is also possible that gravitino is 
the lightest supersymmetric particle, in which case JRJ the upper bound on the reheating 
temperature is 10 11 GeV and thermal leptogenesis is possible. 

We have found interesting correlations between the heavy and light neutrino mass param- 
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eters. For normally-ordered masses, there is an Mn-t dependent upper bound on m Ul , whilst 
for the inverted hierarchy there are both upper and lower bounds on m U3 . Successful lepto- 
genesis with m Ul 0.1 eV is allowed for Mjv x ^ 10 12 GeV. There is also an Mn x -dependent 
lower bound on m ee for normally-ordered light neutrinos, implying its possible testability in 
future experiments. On the other hand, m ee has a preferred value determined by Am 2 tm 
even in the case of the inverted mass hierarchy. It tends to be below O(10 _1 ) eV, making 
the discovery of (3(3q u decay in current experiments improbable. 

There is no correlation between leptogenesis, and the LFV decays and EDMs of charged 
leptons. The branching ratios for r — > ^{e) r y and /i — > may saturate their present lower 
bounds, and the EDMs of electron and muon reach d e ~ io~( 27 ~ 28 ) e cm and ~ 10~ 25 e 
cm in our random samples. There is some correlation between and the LFV r decays 
for normally-ordered light neutrino masses: observation of Br{r — *> cy) would require 
to be an order of magnitude above the lower bound. 



The type of parametrization discussed in |17j and applied in [[18] and in this work is a 
useful tool for studying the large parameter space of the minimal supersymmetric seesaw 
model. In the future, it may help attempts to devise an experimental strategy for determining 
systematically all the 18 parameters in this sector. It would also be interesting to combine 
this approach with models aiming at predictions for some (all) of the seesaw parameters. 
As we have shown in this paper, this type of parameterization can teach us some salutary 
lessons about leptogenesis and its relations to other observables. 
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